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ABSTRACT 



We present ISOPHOT.S, ISOSWS and 8 to 13 mi- 
cron ground based observations of Centaurus A that 
show prominent PAH and silicate features. These 
and other data are used to construct a model for the 
infrared continuum. We find that in a nuclear sized 
aperture (^4 arcsec,^60 pc) the SED is characteris- 
tic of emission from a starburst and AGN torus; in 
larger apertures an additional component of cirrus 
emission is required. Based on our model, the torus 
diameter is estimated to be 3.6 pc and the best fit- 
ting inclination angle of the torus is 45 degrees. This 
result has implications for the detectability of tori in 
low power AGN and in particular for the use of the 
IRAS 60/25 micron flux ratio as an indicator of the 
torus inclination. 
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1. INTRODUCTION 



idence of a merger (see Mirabel et al, these proceed- 
ings), starburst and AGN activity, with HII regions, 
shells, jets, optical filaments and a warped dust lane. 

Evidence of star formation is most apparent in the 
dust lane that intersects the host galaxy NGC5128. 
Marston and Dickens (1988) found that their 12 mi- 
cron IRAS DSD observations followed the Ha emis- 
sion, along the dust lane, and hence the regions of 
star formation. They modelled their 12, 25, 60 and 
100 micron observations as a cirrus spectrum of small 
and large grains heated by the interstellar radia- 
tion with two grain temperatures - hot (240K) small 
grains and cooler (30K) large grains. Evidence of 
starburst activity comes from enhanced far-IR and 
sub-mm nuclear emission that shows a structure off- 
set from the dust lane (e.g. Hawarden et al, 1993) and 
from the high star formation rate, which is typical of 
a starburst galaxy (Eckart et al, 1990). Evidence 
of AGN activity comes from the radio jets/lobes 
(e.g. Clarke, Burns and Norman, 1992) and variable 
x-ray emission (e.g. Morini, Anselmo and Molteni, 
1989) whilst evidence of a Seyfert-like dusty torus 
comes from the high optical depth to the nucleus (e.g. 
Blanco, Ward and Wright, 1990). 



Rowan- Robinson and Crawford (1989) attribute the 
infrared (IR) spectra of galaxies to a mixture of up 
to 3 components: i) general disc emission from grains 
heated by the interstellar radiation field (cirrus), ii) 
a Seyfert component peaking in the mid-infrared and 
iii) a starburst component peaking at about 60 mi- 
crons. Essentially all of these sources arise from the 
thermal reprocessing of ultra-violet and other high 
frequency photons by the dust within these objects. 
The dust grains re-radiate the absorbed energy in 
the IR, with the resultant spectrum dependent on 
the distribution of dust grain temperatures. For the 
starburst and Seyfert components, the clouds of dust 
are optically thick, even to IR photons, so radiative 
transfer effects are important. 

Centaurus A, the famous southern radio (FRI) coun- 
terpart to NGC5128, identified by Bolton, Stanley 
and Slee (1949), at a distance of approximately 3.1 
Mpc (Tonry and Schechter, 1990) is the closest active 
galaxy to us. It is a multi-faceted object, showing ev- 



2. OBSERVATIONS AND DATA 



The observations used to constrain the models are an 
array of published and unpublished data, taken from 
a variety of sources and grouped into 3 aperture sizes: 
a small nuclear aperture of ~4 arcsec, an intermediate 
aperture of ~20 arcsec and a large aperture of "^90 
arcsec. 

The ISOPHOT_S observations (intermediate aper- 
ture) have not been presented before, see figure 1. 
The observations have been reduced with PIA with- 
out additional modifications. The zodiacal light con- 
tributions have been removed using a chopped mea- 
surement. The ISOPHOT_S spectrum consists of two 
parts - PHT_SL (6 to 12 microns) and PHT.SS (2.5 
to 5.5 microns). Clearly seen in the spectrum are 
the deep 9.7 micron silicate absorption feature and 
PAH emission features, however importantly, PAH 
emission at 3.3 microns is not clearly detected. 
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Figure 1. ISOPHOT.S observations. The flux 
units are ergs~^ cm~^ . The wavelength units are 
LOG (microns). 

The ISOSWS observations (intermediate aperture) 
presented here are the continuum points. The emis- 
sion hne spectrum, details of the observations and 
data reduction wiU be presented in Sturm et al (in 
preparation) . The continuum points were determined 
by a second order polynomial fit to the line free re- 
gions of the data. The error bars represent the un- 
certainty in flux calibration - the RMS noise level is 
insignificant by comparison. No correction has been 
made for the different aperture sizes at the different 
wavelengths. 

Due to the problem of starlight contamination, and 
some uncertainty in the flux levels at the short wave- 
length end of these ISO spectra, the data shortward 
of 6 microns are not used in the model fitting. 

The UCL spectrometer observations (nuclear aper- 
ture) have not been presented before, see figure 2. 
The silicate feature at 9.7 microns is clearly detected. 
In addition there appears to be a component of 11.3 
micron PAH and 12.8 micron [Nell] emission. 

Additional data was obtained from the literature. 
The near-IR points (nuclear aperture) were taken 
from Packham et al (1996) and have had the stellar 
contribution removed. The sub-mm points (interme- 
diate and large aperture) were taken from Hawarden 
et al (1993) and have been corrected for non-thermal 
contributions. The 12, 25, 60 and 100 micron large 
aperture points were taken from Marston and Dick- 
ens (1988) and are IRAS DSD observations. 



3. THE MODEL COMPONENTS 



The proposed model for Centaurus A is a combina- 
tion of a number of components. The first component 
is the AGN dusty torus, based on the NGC1068 ta- 
pered disc infrared model of Efstathiou, Hough and 
Young (1995) which assumed a quasar-like central 
continuum source illuminating an optically and ge- 



Figure 2. UCL spectrometer observations. The units 
are the same as in figure 1. 



ometrically thick dusty torus. The main model pa- 
rameters are: opening half-angle of the toroidal cone 
(6>), ratio of the inner and outer radii (rl/r2), ratio of 
the height to outer radius (h/r2), equatorial optical 
depth {cruv)^ dust sublimation temperature (Tl) and 
radial dependence of the density distribution (r~^). 
The second and third components are the starburst 
and cirrus models of Efstathiou, Rowan-Robinson 
and Siebenmorgan (2000). The cirrus spectrum is 
assumed to arise from 12,000 K stars centrally illumi- 
nating a diffuse optically and geometrically thin dust 
shell. The starburst is modelled as an ensemble of 
giant molecular clouds centrally illuminated by hot 
stars. The stellar population is modelled in terms 
of the population synthesis models of Bruzual and 
Chariot (1993) and the giant molecular clouds evolve 
with time due to the expansion of the HII regions. 
Additional information on the model components is 
given in Alexander et al, 1999. 

To take into account of dust lane absorption, radi- 
ation from the torus and starburst are visually ex- 
tincted by 10 mags, essentially the same as the value 
adopted by Packham et al (1996). The cirrus is not 
extincted since it is itself associated with the dust 
lane. The extinction model used assumes the grain 
model and interstellar extinction curve of Rowan- 
Robinson (1992). 



4. THE CENTAURUS A MODEL 



For small apertures, and at IR wavelengths below 
12 microns, it was assumed that the SED would be 
dominated by radiation from the torus, seen through 
the dust lane. Therefore the torus model is fitted to 
this wavelength range, which includes the 9.7 micron 
silicate feature. In fitting to larger apertures, a con- 
tribution from the smaller apertures is included and 
any additional components required. 

The torus parameters were set to be the same as those 
in the Efstathiou, Hough and Young (1995) model 
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for NGC1068 {0=30 degrees, rl/r2=0.01, h/r2=0.1, 
e^^^=l,200, Tl=950 and f3=-l) and the SED was 
calculated for different values of inclination (the an- 
gle between the polar axis and the line of sight). No 
value of inclination gave a simultaneously good fit to 
both the near-IR continuum and the silicate feature, 
with the former best fit by 40 degrees and the latter 
by 70 degrees. The observed silicate feature shows a 
feature at 11.3 microns which is probably PAH emis- 
sion, therefore a component of starburst was added. 
This provided an excellent fit to the silicate feature, 
with roughly equal components of torus and starburst 
for wavelengths less than 60 microns, with the torus 
at an inclination of 45 degrees, see figure 3. 

For larger apertures a component of cirrus emission 
was required. For completeness, the contribution this 
would make in the nuclear aperture is shown in figure 
3 (the cirrus component was determined assuming 
constant surface brightness, as implied from the fits 
to the larger apertures). 
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Figure 3. The model fit to the nuclear aperture data 
points. At the 100 micron wavelength: the lowest 
intensity model component is the cirrus, the highest 
intensity component is the starburst and the compo- 
nent in between is the torus. The units are the same 
as in figure 1. 

In the intermediate aperture, the torus component is 
kept fixed and an additional contribution from the 
starburst is added, producing a reasonable fit to the 
intermediate aperture data points for the silicate fea- 
ture and IR continuum below 25 microns. The PAH 
features and the continuum at longer wavelengths 
showed that an additional component was required. 
Observations (e.g. Marston and Dickens, 1988) sug- 
gest that this additional component is cirrus emission 
from stars within the dust lane. By slightly reduc- 
ing the contribution from the starburst and adding a 
component of cirrus, an excellent fit to all the data 
points is obtained, see figure 4. 

The observed size of the starburst component (e.g. 
Hawarden et al, 1993) limits the contribution it 
should have in the large aperture. Whilst there 
should be an additional contribution from the star- 
burst, the cirrus emission, which is associated with 
the dust lane, and so is widespread, should contribute 
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Figure 4- The model fit to the intermediate aperture 
data points. The model components and units are the 
same as in figure 3. 



significantly more. The best fit is obtained with a rel- 
atively small increase in the starburst contribution 
and a significant increase in the cirrus contribution, 
see figure 5. 
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Figure 5. The model fit to the large aperture data 
points. The model components and units are the same 
as in figure 3. 



5. DISCUSSION 



The complete model includes contributions from 
torus, starburst and cirrus in all apertures. The 
model fit in the nuclear aperture is quantitatively 
similar to that found by Laurent et al (these proceed- 
ings), with ISOCAM CVF observations and AGN 
and starburst observational templates. 

The opening half-angle of the torus cone is very sim- 
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ilar to that of the 'cones' in the recently pubhshed 
HST nuclear Pace image (Schreier et al, 1998), and 
if it is assumed that these observations are of ionisa- 
tion cones, rather that the suggested warped accre- 
tion disc, the cone opening half-angle is constrained 
between 30 and 40 degrees. Extinction to the 2.2 mi- 
cron region can be estimated by comparing the flux 
for the face-on torus to the 45 degree inclined torus 
and is 2.6 mags, equating to a visual extinction of 23 
mags. This degree of extinction is consistent with the 
analysis of Meadows and Allen (1992) and suggests 
that the near-IR emission region is considerably less 
extincted than the x-ray region (70 mags, Blanco, 
Ward and Wright, 1990). 

Based on this model, the torus diameter is calculated 
to be 3.6 pc, assuming a distance of 3.1 Mpc, placing 
it easily within the nuclear aperture. This is signif- 
icantly smaller than that actually observed in other 
galaxies (e.g. "^1 kpc for MG0414+0534 (Oya et al, 
these proceedings), '^200pc for NGC1068 (Young et al, 
1996), -^100 pc for NGC4261 (Jaffe et al, 1993) and 
-^50 pc for NGC4151 (Mundel et al, 1995)). How- 
ever, although of a small diameter, this model does 
not support the even smaller diameter tori of the Pier 
and Krolik (1992) model. 

Possible observational support for the small torus size 
estimate comes from the 2.3 GHz radio observations 
of the nuclear sub-arcsecond scale jet and counter jet 
(Jones et al, 1996). These images show the core to be 
completely absorbed between the jet and counter jet, 
attributed to a gaseous disc or torus of 0.4 to 0.8 pc. 
The model torus, with an inclination of 45 degrees 
would present a 1.4 pc absorption band; very similar 
to the size of the radio core absorption. 

Within the context of unified theories the inclination 
of the torus is a crucial parameter in determining 
whether an AGN is a Type 1 or Type 2. Various stud- 
ies have suggested that the IRAS 60/25 micron flux 
ratio is an indicator of torus inclination with high or 
low ratios implying edge-on to face-on torus inclina- 
tions. The combination of all the model components 
gives a 60/25 micron flux ratio of 9.4, although the 
actual torus 60/25 micron flux ratio is only 0.7. The 
torus inclination implied from this study is identical 
to that found for NGC1068 (Efstathiou, Hough and 
Young, 1995) and yet the IRAS 60/25 micron flux 
ratio for NGC1068 is only 2.1. If the torus was as 
powerful as that in NGC1068, or a smaller 'nuclear' 
sized aperture was used, the 60/25 micron flux ratio 
would be more indicative of the torus inclination but 
within the large IRAS sized aperture, the torus is 
dominated by other more powerful IR emitting com- 
ponents. 

A comparison of the mean IRAS 60/25 micron flux 
density ratios of the extended 12 micron sample of 
galaxies, which is considered statistically complete on 
12 micron flux, gives a mean of 4. lib 2. 9 for the Seyfert 
Is and 4.5=b3.0 for the Seyfert 2s (Andy Thean, pri- 
vate communication). This suggests that this colour 
ratio cannot distinguish between Type 1 and Type 
2 AGN and so cannot be a reliable indicator of the 
torus inclination. 

The IRAS 60/25 micron flux density ratio is more 
likely to provide a ratio of the AGN to star forma- 
tion components. Evidence for this comes from the 
60 micron to 6 cm flux correlation (Wilson, 1988) 



which works only if the radio jet/lobe AGN emis- 
sion is removed (Peter Bartel, private communica- 
tion) , thereby strongly suggesting that the 60 micron 
flux is coming from star formation and starburst re- 
gions. Additional evidence comes from the ISOLWS 
observations of NGC1068 (Spinoglio et al, these pro- 
ceedings) where the far-IR continuum is almost in- 
distinguishable to that of the archetypal starburst 
galaxy M82. 
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